Differential rotation enhanced dissipation of tides in the PSR J0045-7319 Binary 



Pawan Kumar * and Eliot J. Quataert 
Department of Physics 
Massachusetts Institute of Technology, Cambridge, MA 02139 

Abstract 

Recent observations of PSR J0045-7319, a radio pulsar in a close eccentric orbit with 
a massive B-star companion, indicate that the system's orbital period is decreasing on a 
timescale of ~ 5 x 10 5 years, which is much shorter than the timescale of ~ 10 9 years given 
by the standard theory of tidal dissipation in radiative stars. Observations also provide 
strong evidence that the B-star is rotating rapidly, perhaps at nearly its break up speed. 
We show that the dissipation of the dynamical tide in a star rotating in the same direction 
as the orbital motion of its companion (prograde rotation) with a speed greater than the 
orbital angular speed of the star at periastron results in an increase in the orbital period 
of the binary system with time. Thus, since the observed time derivative of the orbital 
period is large and negative, the B-star in the PSR J0045-7319 binary must have retrograde 
rotation if tidal effects are to account for the orbital decay. We also show that the time 
scale for the synchronization of the B-star's spin with the orbital angular speed of the star 
at periastron is comparable to the orbital evolution time. From the work of Goldreich 
and Nicholson (1989) we therefore expect that the B-star should be rotating differentially, 
with the outer layers rotating more slowly than the interior. We show that the dissipation 
of the dynamical tide in such a differentially rotating B-star is enhanced by almost three 
orders of magnitude leading to an orbital evolution time for the PSR J0045-7319 Binary 
that is consistent with the observations. 
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1. Introduction 



An unusual binary system was recently discovered in the Small Magellanic Cloud 
(SMC); Johnston et al. 1994, Kaspi et al. 1994. One of the members of this binary is a 
pulsar, PSR J0045-7319, with a spin period of 0.93 s; its companion is a main sequence B- 
star of mass ~ 8.8M©. The orbital period of this binary system is 51 days, the eccentricity 
is 0.808, and the separation between the stars at periastron is about 4 times the B-star 
radius. Kaspi et al. (1996) report that the orbital period for this system is evolving on a 
time scale of ~ 5 x 10 5 years, whereas the standard theory for tidal dissipation predicts an 
evolution time of order 10 9 years. Observations also indicate that the spin of the B-star is 
not perpendicular to the orbital plane and that the component of the spin normal to the 
orbital plane is close to the break up rotation speed (Kaspi et al. 1996, Lai et al. 1995). 

It has been suggested (Lai 1996) that the rapid evolution of the orbital period of 
the SMC binary can be accounted for by the radiative dissipation of the dynamical tide, 
provided that the B-star companion to PSR J0045-7319 has a rapid retrograde rotation 
rate. This is an interesting idea and the observations do suggest that the B-star is rapidly 
rotating, though the sign of the spin is unknown. We show, however, that as long as the B- 
star is rotating rigidly, the dissipation of the dynamical tide is too inefficient to explain the 
observed evolution of the orbital period. When the star has some differential rotation, so 
that the frequency of the tidally excited waves in the rest frame of the fluid near the surface 
of the star is smaller than in the deep interior, the wave dissipation increases dramatically, 
readily explaining the observations. A more detailed and rigorous presentation of this 
analysis will be given in a future paper (Kumar & Quataert 1996). 

2. The dynamical tide in a rigidly rotating star 

Waves excited by the tidal forcing have frequencies ~ 2Q p as seen by an inertial 
observer, where fi p is the angular speed of the star at periastron. If the angular velocity 
of the B-star perpendicular to the orbital plane is Q*, then the wave frequency in the local 
rest frame of the fluid is ujt* ~ 2[fi p — O*]. The angular momentum luminosity in the 
dynamical tide (the angular momentum transported by waves across a sphere of radius r 
per unit of time) is 

Ltide = 2pr 2 uu^t{r)v g , (1) 

where £f is the square of the transverse displacement (integrated over solid angle) associ- 
ated with the dynamical tide, v g ~ rui^/ (6 X / 2 A^) is the group speed of quadrupole gravity 
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waves in the limit of small wave frequencies {oj t * <C N), and N is the Brunt-Vaisala fre- 
quency. We note that the radial displacement amplitude of the dynamical tide is much 
smaller than its transverse displacement and therefore has been neglected in the above 
equation. It can be shown (Zahn 1975) that is approximately equal to the transverse 
displacement associated with the equilibrium tide, i[ eq \ about a wavelength away from the 
convective core. Furthermore, it is easy to show that the squared transverse displacement 
of the equilibrium tide integrated over all angles is 
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(2) 



1.4M is the mass of the neutron star, 
and g is the gravitational acceleration of the primary star. The total angular momentum in 
the dynamical tide is given by L t ide ~ T per jL ti d e , where T peri ps 7r(l — e) 2 / (0 o Vl — e 2 ) ~ 2 
days is the duration of the periastron encounter and e = 0.808 is the orbital eccentricity. 
The rate of change of the orbital energy due to the dissipation of the dynamical tide is 
E or b = —^E tide , where V is the dissipation rate for the dynamical tide (calculated below) 
and Eude is the energy in the dynamical tide, which can be shown to be equal to L^ e fi p , 
i.e., 
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where M t ss 10. 2M© is the total mass of the binary system, fi G = (GMt/a 3 ) 1 / 2 = 1.42 x 
10~ 6 rad s -1 is the orbital frequency and a 126R & is the semi-major axis of the orbit. 
Note that all quantities in the above equation that depend on r are evaluated about one 
wavelength away from the inner turning point of the wave. We take the radius of the B-star 
to be 6.0.R©; the radius of the convective core of the B-star, where the gravity waves become 
evanescent, is about 0.23 .R*, the mass in the convective core is ^ 3M and the density 
and Brunt-Vaisala frequency outside the core are 2 g cm -3 and 100 /uHz, respectively. 
Substituting these into equation (3) we find that the total energy in the dynamical tide for 
a non-rotating star, so that uj t * ~ 20 p ps 7.2[iRz, is ps 10 40 ergs. For retrograde rotation 
of the B-star at a frequency of fi*/27r = -8^Hz, or fi* = fi^i^/GM*) 1 / 2 ps -0.4, the 
frequency of the tidally excited wave is u t * ~ 23/uHz and equation (3) implies that the 
energy in the tide is a factor of about 100 larger compared to the non-rotating star. This 
point was made by Lai (1996) who suggested that the enhanced energy in the dynamical 
tide can explain the observed evolution of the SMC system. Unfortunately, as we point 
out below, this is not tenable because the damping time of a wave of frequency ~ 23/uHz 
is also larger by a factor of about 100 than a wave of frequency ps 7.2(j,}Iz. 
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For prograde rotation the tidal frequency and energy decrease with increasing rotation 
rate so long as < O p (see eq. [3]). The tidal energy does not actually vanish when 
fi* = O p because axisymmetric waves are still excited. For fi* > fi p the tidal torque causes 
the angular momentum of the stellar spin to decrease (which is taken up by the orbit) and 
as a result the orbital energy and the orbital period increase with time. Equation (3) shows 
that Ende < 0, and so E or b > 0, when fi* > Q p , which may seem paradoxical. We note 
that the tidal energy as seen in the rotating frame of the star is a positive definite quantity; 
the tidal energy as seen from an inertial frame (given in eq. [3]), however, also includes 
a contribution from the work done by the tidal force on the rotating fluid in the star, 
which can either increase or decrease the star's rotational kinetic energy. The negative 
tidal energy for fi* > fi p corresponds to more energy being taken out of the rotation than 
is deposited in the dynamical tide. We emphasize that the tidal energy is positive so long 
as < O p , which obviously includes the case of arbitrary retrograde rotation, and is 
negative when Q* > O p . Thus, if the rotation rate of the B-star is greater than O p , the 
observed period decrease of the PSR J0045-7319 binary implies that the direction of the 
rotation of the B-star must be retrograde. We also note that the energy in the tide for 
O* O p is too small by about 2 orders of magnitude to explain the observed rate of period 
evolution. 

A detailed calculation of the damping of gravity waves is given in Kumar & Quataert 
(1996). Here we present a crude estimate to point out that the wave damping increases 
very rapidly with decreasing wave frequency. 

The local damping rate of a wave, 7(r), is approximately equal to the inverse of the 
time it takes photons to diffuse across a wavelength times the ratio of the photon energy 
density to the total thermal energy density of the plasma. This can be written in the 
following convenient form, 

7(r)~Ffc r 2 /( W ), (4) 

where F is the radiative flux of the star and k r ~ 1(1 +1)/ (ru t *) is the radial wave 
number for gravity waves (£ — 2 for quadrupole waves) . We see from the above expression 
that the wave damping increases with decreasing density in the star and thus most of 
the contribution to the wave damping comes from near the upper turning point of the 
wave which is located at a depth of ~ R^uj^/(6g) beneath the stellar surface. Integrating 
the local dissipation rate, weighted by the wave's energy density, yields the global wave 
damping rate. Since the energy density in the wave is inversely proportional to its group 
speed, we find that the global damping rate for the wave, in a star of polytropic index 
2 (a value appropriate for the outer envelope of the Yale group's B-star model we are 
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using), scales as cu t ~ 7 . This scaling only applies to waves with more than about 4 radial 
nodes, which corresponds to u) t * ^ 15/uHz for the SMC B-star. Modes with fewer nodes 
or higher frequencies have wavelengths of order the stellar radius and their dissipation is 
not as concentrated near the outer turning point as it is for lower frequency modes. Thus 
the global damping rate of high frequency g-modes has a much weaker dependence on 
mode frequency. Using equation (4) and the parameters for the SMC B-star given above 
we estimate that the damping time for a low order quadrupole gravity wave is ~ 3000 
years, which is consistent with the results of a number of different nonadiabatic eigenvalue 
calculations (Kumar & Quataert 1996; Saio & Cox 1980). 

The frequency of the dynamical tide in the non-rotating star is ps 7.2/jH.z, while for 
rapid retrograde rotation (fi* = —0.4) the frequency is w 23/uHz. Using the above scaling 
for the damping of the dynamical tide, we find that the ratio of the damping times in 
these two cases is 1/200. Since the rate of dissipation of the orbital energy is equal to 
the tidal energy times its damping rate, we see that even though the energy in the tide 
for rapid retrograde rotation is a factor of ~ 100 larger than for no rotation, the rate of 
dissipation of the orbital energy is very nearly the same for fi* = —0.4 and O* = 0. Using 
the dissipation time for the dynamical tide in the retrograde rotating star of ~ 3000 years, 
the energy in the dynamical tide of ~ 10 42 ergs, and the orbital energy for the SMC binary, 
E orb ^ 2 x 10 47 ergs, we find that the orbital period evolution time for the SMC Binary is 
^ 10 9 years, or roughly three orders of magnitude larger than the observed value. 

3. Differential rotation enhanced dissipation of tides 

It should be clear from this discussion that one way to get more rapid evolution of the 
orbit is to find a more efficient way of dissipating the energy in the dynamical tide. This 
would occur naturally, considering the sensitivity of wave dissipation on frequency, if the 
frequency of the tidally excited waves were smaller near the surface of the star, where the 
dissipation occurs. Waves with frequencies ~ 20 p , as seen in an inertial reference frame, 
are excited near the interface of the convective core and the radiative exterior where the 
wavelength of the wave is the largest. Thus, in order to get a smaller wave frequency in 
the rest frame of the fluid near the surface of the star, we need the surface to be rotating 
at a frequency close to Vt p ~ 3.6/uHz. We show below that this is in fact expected for the 
B-star of the SMC binary. 

For nearly circular orbits, the rotation of stars in close binary systems tends to syn- 
chronize with the orbital motion, that is, the rotation frequency of the star approaches 
the orbital frequency. For highly eccentric orbits, the tidal force is only appreciable near 
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periastron and so the rotation frequency of the star approaches Q p . This is called pseudo- 
synchronization. 

We estimate the relative timescales for spin pseudo-synchronization and orbital cir- 
cularization for the SMC system to determine if the B-star should be rotating pseudo- 
synchronously. The pseudo-synchronization timescale due to the transfer of orbital an- 
gular momentum to the star is L*/L* = — (L*/L orb )(L orb /L orb ), where L orb /L* ps 4fi~ 1 
for the SMC Binary and L orb E orb /Q p (see above). Since L orb Q /E orb w 1, we find 
that L*/L* ps 5£l*E orb / 'E orb for the SMC Binary system. It follows from the general 
equations for tidal evolution in the weak friction limit that a similar result holds for the 
equilibrium tide (Hut 1981). We conclude that so long as the B-star in the SMC Binary 
is rotating near break up, both the equilibrium and the dynamical tides lead to timescales 
for spin pseudo-synchronization which are comparable to the orbital circularization time. 
We therefore do not expect that the B-star's interior has spun down appreciably, which is 
clearly consistent with the observations. 

We do expect, however, that the surface of the B-star will have slowed down, leading 
to appreciable differential rotation in the B-star. The physical reason for this follows from 
the seminal work of Goldreich and Nicholson (1989). Waves generated by tidal forces 
deposit their angular momentum at the place in the star where they are dissipated. Since 
the dissipation rate is largest near the surface of the star, and the moment of inertia of the 
surface region is a small fraction of the star's total moment of inertia, the outer region of 
the star tends to be pseudo-synchronized on a much shorter time scale than the interior. 
For the retrograde rotating B-star with O* ps —0.3, the frequency of the dynamical tide 
is ~ 17/uHz; the outer turning point of this wave is at r ps 0.9/?* and so the rotation 
of the B-star at r ^ 0.9.R* should be pseudo-synchronous. Those parts of the star lying 
above or significantly below this radius, however, experience little net tidal torque and thus 
can continue to rotate at the primordial rotation rate of the star provided that magnetic 
stresses and shear instabilities do not efficiently transport angular momentum to/from the 
pseudo-synchronous layer. The optical linewidth measurement of Bell et al. (1995) gives 
a projected rotation velocity of the B-star of 113±10 km s _1 . If correct, this suggests that 
the surface of the star is perhaps not pseudo-synchronized, supporting our simple model. 

The frequency of the tidally excited gravity wave in the local rest frame of the 
star a distance r from the center is ut*(r) ~ 2[Q P — 0*(r)] and its wavenumber is ~ 
Q 1 / 2 N(r)/(ruJt*) 1 where 0*(r) is the local angular rotational velocity of the star normal to 
the orbital plane. Thus, as the wave approaches the pseudo-synchronously rotating surface 
of the star its frequency and wavelength approach zero and the energy it carries is entirely 
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dissipated. The energy in the dynamical tide for retrograde stellar rotation at a frequency 
of ~ — 6/uHz (fi* Rs —0.3) is ~ 10 41 ergs (see eq. [3]). With a change in the orbital energy 
per orbit of AE orb —E tide xs — 10 41 ergs, and the orbital energy of 2 x 10 47 ergs, we 
find that P or b/P rb = 1.5E or t/E or t ~ 4.2 x 10 5 years, which in good agreement with the 
observations. It is easy to show that the rate of change of orbital eccentricity (e) is a factor 
of about 7 smaller than P or b/P or b- Thus the theoretically expected value for e is a factor 
of ~ 20 smaller than the current observational limit (Kaspi et al. 1996). 

We note that prograde rotation of the B-star at O* ~ 0.4 gives an orbital evolution 
time of comparable magnitude, though with the opposite sign. This emphasizes that 
retrograde rotation is not required to explain the observed short timescale for orbital 
evolution, but rather to get the correct sign for the time derivative of the orbital period. 
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